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ABSTRACT 

bJQi We use the progenitor of SN 2012aw to illustrate the consequences of modeling circum- 

stellar dust using Galactic (interstellar) extinction laws that (1) ignore dust emission in the 
near-IR and beyond; (2) average over dust compositions, and (3) mis-characterize the opti- 
cal/UV absorption by assuming that scattered photons are lost to the observer. The primary 
consequences for the progenitor of SN 2012aw are that both the luminosity and the absorption 
are significantly over-estimated. In particular, the stellar luminosity is most likely in the range 
10 4 8 < L*/L & < 10 5 and the star was not extremely massive for a Type IIP progenitor, with 
M* < 15M Q . Given the properties of the circumstellar dust and the early X-ray/radio detec- 
tions of SN 2012aw, the star was probably obscured by an on-going wind with M ~ 10~ 5 5 to 
lO~ 5 O M /year at the time of the explosion, roughly consistent with the expected mass loss rates 
for a star of its temperature (T* ~ 3600^2(jo K) and luminosity. In the spirit of Galactic extinction 
laws, we supply simple interpolation formulas for circumstellar extinction by dusty graphitic 
and silicate shells as a function of wavelength (A > 0.3fim) and total (absorption plus scatter- 
ing) V-band optical depth (ry < 20). These do not include the contributions of dust emission, 
but provide a simple, physical alternative to incorrectly using interstellar extinction laws. 

Subject headings: stars: evolution - supergiants - supernovae:general 



1. Introduction 

A key component to understanding supernovae (SNe) is the mapping between the explosions and their 
progenitor stars. Slow, steady progress is being mad e, and it is n ow well established that Type IIP SN are 



associated with red supergiants (see the review by ISmartt 2009 ). There is a puzzle, however, in that the 



observed upper limit of ~ 16M Q on the masses Type IIP SN progenitors appears to be signifi cantly lower 



than t he maximum masses of ~ 25M for stars expected to explode wh ile still red supergiants (ISmartt et al. 



2009), part of a more general absence of massive SN progenitors (see iKochanek et al] |2008). Since these 
"missing" progenitors should be more luminous than those which are being discovered, they must either 
be hidden, evolve differently than expected, or fail to explode. For example, in the rotating models of 
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Ekstrom et al.1 (120121) . the upper mass limit for red s upergiants is lower, with IQ M^ stars being blue r ather 
than red at the onset of carbon burning. Alternatively, lO' Connor & Ota (1201 ll) and lUgliano et all (120121 ) find 
that the progenitor masses of 20-25M Q corresponding to the upper mass range for red supergiants are more 
prone to failed explosions and prompt black hole formation, and such events would ha ve to be found by 
searching for stars disappearing rather than explosions appearing (IKochanek et al.l l2008). Binary evolution 
can also alter the distribution of final states at the tim e of explosion given the high probability for mass 
transfer as stars expand to become red supergiants (e.g. ISana et al.ll2012r) . 



If the discrepancy is not explained by the physics of stellar evoluti on or explosion, then the most likely 
remaining explanation is the effect of circumstellar dust. For example, IWalmswell & Eldridgd (120121 ) note 
that more massive and luminous red supergiants have stronger winds which can form dust and partly obscure 
the star. This then biases (in particular) the upper mass limits associated with failed searches for progenitor 
stars, since the lum inos ity limits must be corr ected for an unknown amount of circumstellar dust extinction. 
Fraser et all (120121) and I Van Dyk et all (120121) recently analyzed the progenitor of S N 2012aw, finding that it 
was both relatively high mass (15-2OM for lvan Dyk et al.l(|2012l) and 14-26M for lFraser et aD(|2012lv ) and 
the most heavily obscured of any SN progenitor other than the comp letely obscured (and debated) SN 2008S 
class (see IPrieto et al.ll2008l iThompson et al. 2009, Koch anekll201 ll). Since most of the exti nction vanished 
after the SN, the dust must have been circumstellar (IFraser et al.ll2012llVan Dyk et al.ll2012l) . 



Li ke most studies of c ircu mstellar dust in superno vae or supernovae progenitors, IWalmswell & Eldridge 



<l2012l) . 



Fraser et al. 



(|2012|) and lVan Dyk et al.1 (|2012) treat circumstellar dust as if it is a foreground screen 
that can be quanti tatively modeled us ing Galactic interstellar extinction curve models parametrized by the 
value of Ry (e.g. ICardelli et al.lll989i) . Usually only studies that are self-consistently calculating emission 



by the dust correctly model the absorption by ci rcu mstellar dust (e.g. stu dies of the SN 2008S class of 
transients by IWesson et al.ll201Ql lKochanekll201 ll or lSzczygiel et al.ll2012r) . But three well-known effects 



mean that it is never appropriate to make this approximation unless the optical depth is negligible compared 
to the required precision of the analysis. First, emission from circumstellar dust can be important in the 
near-IR if the star is presently forming dust at temperatures of 1000-2000 K. This matters most for hotter 
stars with less intrinsic near-IR emission than the red supergiants we consider here. Second, interstellar dust 
has the average composition of dust from all sources, while individual stars have the dust associated with 
their particular che mistry. This is relevan t here because ~ 20M Q stars generally produce silicate rather than 
graphitic dusts (e.g. IVerhoelst et al.l l2009). 



The third, and least appreciated point, is the different role of scattered pho tons in inters t ellar a nd cir- 
cumstellar extinction. This issue has been discussed in detail by IWangl ((2005) and iGooban (120081) in the 
context of anomalously low estimates of Ry for Type la SN, but correctly treating the problem has not 
become a matter of practice. For a foreground screen, scattered light forms a very diffuse, extended halo 
around the source which is not included in the estimate of the source flux. This halo can sometimes be seen 



as the extended dust echoes of transient sources (e.g. ISugerman & Crottsll2002l) . If the circumstellar dust is 
unresolved, however, the scattered light is simply included in the total flux. Thus, for interstellar extinction 
you observe only the direct, unabsorbed and unscattered emission, while for circumstellar extinction you 
observe both the direct unabsorbed and the scattered emission. Most of the observed optical emission from 
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Fig. 1. — DUSTY model of a 10 4 K black body surrounded by ry = 3 of cold silicate dust in a shell 
with Rout/Rin = 2 and p oc 1/r 2 in the shell. The dust converts the input spectral energy distribution (SED) 
to the output SED, where the output SED is comprised of contributions from photons that escape without 
scattering (direct light) and photons that escape after being scattered (scattered light). In an image, the direct 
emission is a point source and the scattered emission is a halo with a radius comparable to the inner radius 
of the shell. If the dusty shell is large enough to be resolved, only the direct emission is measured as coming 
from the source. We have made the dust cold enough to have no contributions from dust emission over this 
wavelength range. 
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Fig. 2. — Effective extinction laws R\ if the dust is a foreground screen where we only measure the 
unscattered direct light. The solid curves show the results for ry = 3 of graphitic dust ( top), silicate dust 
(bottom) and a 50:50 mix (middle). The dashed curves show the ICardelli et all (119891) extinction laws 
(CCM) with R v = 4 (top), 3.1 (middle) and 2.0 (bottom). Normal Galactic dust (R v = 3.1) is relatively well 
modeled by the 50:50 mix. 
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Fig. 3. — Effective extinction laws R\ for dust in an unresolved shell around the star where we measure 
both the direct and scattered light. The solid curves show the results for -ry = 3 of graphi tic dust (top), 
silica te dust (bottom) and a 50:50 mix (middle). The dashed curves show the best fit CCM Jcardelli et al. 
1989) extinction laws with Ry = 2.4 (top, graphitic), Ry = 2.1 (middle, 50:50 mix) and Ry = 1.6 (bottom, 
silicate). The fits are never very good, particularly at high o ptical depth, and the best fit parameters vary 
somewhat with the optical depth and wavelength. iGoobarl (120081) obtained similar results. The peak at 
long wavelengths (small A -1 ) is a silicate emission feature which DUSTY includes in the scattered light 
contribution to the emerging spectrum. 
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a moderately self-obscured star is scattered light. 

Fig.[T]shows an example generated by DUSTY ( Ivezic & Elitzurill997 . Ivezic et all 1999 . Elitzur & Ivezic 

2001) using the spectral energy distribution (SED) of a 10 4 K black body surrounded by t v = 3 (scattering 



plus absorption) of cold silicate dust. The dusty material has a density profile of p oc 1 /r 2 in a shell with 
Rout/Rin - 2. Galactic (interstellar) extinction corresponds to putting the dust at such a large radius that we 
no longer include the scattered light in the observed flux of the star, and the extinction law corresponds to 
the wavelength dependent difference between the input spectrum and the directly escaping spectrum. Fig. [2] 
shows this ratio convert ed into an ext i nction law R\ for pure graphitic, pure silicate and a 50:50 mix of 
the two, as compared to lCardelli et all (119891 ) Galactic extinction laws with Ry = 4.0, 3.1 and 2.0. Galactic 
dust m odels typically have a roughly 50:50 mix of graphitic and silicate grains (see the summary in lDraine 
201 lh . and using this mix in the DUSTY models roughly reproduces a typical Ry = 3.1 Galactic extinction 
law. The match is not perfect because there are differences in the assumed size distributions. 

For a star surrounded by an unresolved shell of dust, however, the emission we measure is the sum 
of the direct and scattered light, which is very different from the direct emission alone. Fig. [3] shows this 
case converted into an effective circumstellar extinction law for the same three dust mixtures. For circum- 
stellar dust, all three examples are now significantly below the Ry = 3.1 curve and a given change in color 
implies a sign i ficant ly smaller change in luminosity. If we fit the absorption in these DUSTY models with 
Cardelli et all (|1989l) extinction models, the best fits for the graphitic, silicate and 50:50 dusts are Ry ~ 2.4, 
Ry ~ 1.6 and Ry ~ 2.1, respectively. These are the best fits for ry = 3 from 0.36/um (U band) to 1/xm. The 
best fits vary modestly (ARy ~ 0.1) with ry and the fitted wavelen gth range, become wo rse for higher op- 
extinction curves 



tical depths, and there a re no tru l y goo d fits. Fig. [3] superposes these ICardel 



on the DUSTY models. iGooba r (2008) found similar estimates for the best 



i et al 



il989h 



Cardelli et all (11989b extinction 



curves to model circumstellar dust and also noted that they are not very good approximations. 

The final point to note in Figs.[2]and[3]is that the dust composition is quantitatively important. Interstel- 
lar dust is a mixture of graphitic and silicate dusts from many sources. Individual stars, however, typically 
only produce silicate or graphitic dusts depending on the carbon to oxygen abundance ratio of the stellar 
atmosphere. To zeroth order, all available carbon and oxygen bond to make CO. Then, a carbon poor (rich) 
atmosphere has excess oxygen (carbon) to make silicate (graph itic) dusts. The atmospheres of ~ 20M Q red 
supergiants usually form silicate dusts (e.g. IVerhoelst et al. I l2009h because they neve r undergo the dredge u p 
phases that enrich the atmospheres of the lower mass AGB stars with carbon (e.g. Ilben & Renzinilll983h . 
Assuming a mixed interstellar composition will generally overestimate the absorption associated with a 
given change in color for silicate dusts and underestimate it for graphitic dusts. 

In this paper we consider the effects of dust on the progenitor of SN 2012aw in more detail. First, in 
$2] we use archival Spitzer IRAC data to set stringent limits on the 3.6, 4.5, 5.8 and 8.0/im fluxes of the 
progenitor, and then model the spectral energy distribution (SED) of the progenitor with DUSTY in order 
to correctly model absorption, scattering and emission from a dust enshrouded star. We highli ght where any 



extinc tion law derived for dust screens at great distances from the source such as the standard iCardelli et al. 
(119891) extinction laws can create problems. Finally, in $3] we summarize the results and provide simple 
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interpolation formulas for extinction by graphitic and silicate circumstellar dust as a function of wavelength 
and optical depth. 



2. Mid-IR Flux Limits, X-ray Fluxes and Models 



As noted by Fraser et al.l (|2012l ). there are no obvious sources in the archival SINGS (|Kennicutt et al 



2003) IRAC images of the region around the progenitor. It is, however, a relatively clean region with the 
emission dominated by unresolve d stars, as sh o wn in Fig. [4] We identified 8 nearb y stars in common between 
the HST I-band image used by Fraser et al.1 (12012!) and IVan Dyk et all (120121) to identify the progenitor 
and the image formed by co-adding the 3.6 and 4.5/im data and used the IRAF geomap/geotran scripts to 
estimate the position of the progenitor shown in Fig.|4] Based on jackknife resampling of the reference stars, 
the nominal uncertainties in the position are small (about 0"05). 

Since we are essentially confusion limited, we estimated flux limits using the grid of 36 apertures 
covering the location of the progenitor and bounded by the nearby brighter stars as shown in Fig. [4] The 
grid spacing was 3"0. We used the IRAlf] APPHOT/PHOT to measure aperture magnitudes with a signal 
radius of 2"4 and a sky annulus from 2"4 to I'.'l. The fluxes were calibrated following the procedures given 
in the Spitzer Data Analysis Cookbook and aperture corrections of 1.213, 1.234, 1.379, and 1.584 for the 
3.6 through S.Ofim bands. Keeping the sky annulus immediately next to the signal aperture minimizes the 
effects of background variations created by the galaxy. We used a 2a outlier rejection procedure in order to 
exclude sources located in the local sky annulus, and correct for the excluded pixels assuming a Gaussian 
background flux distribution. 

We clipped 3 of the grid points in the wings of the brighter sources and the two lowest (negative) 
flux points that lay in an obvious "dark lane", leaving 31 points. Based on one of the two SINGS epochs, 
the variances in the four bands are 10.7, 11.4, 12.8 and 14.2 /^Jy for the 3.6, 4.5, 5.8, and 8.0/um bands 
respectively. This is well above the nominal noise level for the exposure time (0.5, 0.8, 6, and 6 /^Jy) 
because there is stellar emission in the region. The variance after coadding the two SINGS epochs is only 
slightly smaller, as expected for noise dominated by confusion, so we rather conservatively use the variances 
derived from the single epoch as 1<t flux limits. Note that some of the grid points clearly lie on sources in 
the coadded image. These sources are not readily apparent in the four individual sub-images, where their 
detection significance would be roughly 2a. A source of similar flux at the supernova position would be 
equally obvious in Fig. [4] 



For m odels of the spectra l energy distribution, we combined the V and I band estimates from lVan Dyk et al 



(|2012l) and lFraser et al.1 (|2012l) . including the difference in t heir V -band estimates as an additional systematic 
error. We used the J and K 4 estimates from lVan Dyk et al.l (|2012l) . based on their better photometric calibra- 
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Fig. 4. — Combined 3.6 and 4.5/mi images of the region surrounding the progenitor of SN 2012aw. The 
2"0 diameter circle marks the estimated position of the progenitor, although the formal uncertainties in the 
position are far smaller (0"05). The smaller points show the grid of apertures used to estimate the flux limits, 
where missing points show the locations of dropped apertures (see text). A few of the grid points lie near 
fainter sources that are obvious in this co-added image but only marginally detectable (at about 2a) in the 
four individual images combined to make Fig. 01 The image is aligned to the world coordinate system with 
North up and East left. 
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tions. We use a minimum flux uncertainty of 20%, which broadens the V, I and J-band photometric errors, 
given that the data were obtained over a long time p eriod (1994-2009, with the SINGS IRAC data being 
obtained on 22/23 May 2004) and lvan Dvk et all J2012h argue for the detection of variability at V and I. We 
model the star using the MARCS (IGustafsson et al.l l2008) stellar atmosphere models , considering only the 
15M , solar metallicity, 5 km/s turbulent ve locity log 10 g = -0.5 models preferred by I Van Dyk et all (120121) 
and also considered by Fraser et all (|2012h . These are available for effective temperatures of T e = 3300 
to 4000 K in steps of 100 K and then 4250 and 4400 K. For intermediate temperatures we linearly in- 
terpolate between the available models. We averaged the high resolution MARCS models and the filter 
transmission functions onto a coarser wavelength grid for use in the DUSTY models and then generated 
magnitude estimates for each model including the appropriate averages over the filter band passes. We first 
fit the measurements to normalize the luminosity, and then added a contribution from the upper limits as 
Ax 2 = (Fmodel/Flimit) 2 f° r eacn band, where F moi i e i is the estimate from the normalized model and Fu m i t is 
the la limit on the flux. 



Imm 



week and 



Stockdale et al. 



er & Brown (2012) reporte d a Swi ft (Gehr els et all 120041) X-ray detection of the SN in the first 



(|2012f) and lYadav et all (|2012l ) report 20 GHz radio detections in the first month. 
Swift continued to monitor the SN for an extended period, so as an independent probe of the circumstellar 
medium we obtained all the archival Swift XRT data, binning the observations as summarized in Table [T] 
We reprocessed the XRT data using the xrtpipeline tool provided by the Swift team. We binned the 
observations into six roughly logarithmic time intervals and then reprojected each group of observations into 
a single image. We chose the source region to be a circle centered on the SN with a radius of 10 pixels (23 "6) 
and a nearby background region w ithout any s ources, using aperture cor rections based on the analysis of the 



Swift PSF by Moretti et all (120051). Following 



Immler & Brownl (|2012l) we assumed a G alactic foreground 



colum n density of 2.9 x 10 20 (|Dickey & Lockmanlll990h . For our epoch best matching llmmler & Brown 
J2012h we reproduce their count rates. While none of the X-ray detections are of very high significance, 
there appears to be a low level of detectable (3a) and probably time varying X-ray emission for the first 1-2 
months. There are no clear detections at low energies (0.2-0.5 keV) and only very marginal detections at 
high energies (2-10 keV) - the observed counts are completely dominated by the 0.5-2 keV band. 



We use the DUSTY jlvezic & Elitzurlll99l llvezic etal.lll999L Elitzur & Ivezic1l2001|) dust radiation 
transfer models to correctly include dust emission, absorption, scattering an d composition when modeling 
the proge nitor SEP . We u se either graphitic or silicate dust models from Draine & Led (| 19841) and the 
default Mathis et al.1 (1 19771) power-law grain size distribution (dn/da oc a -3 5 for 0.005/mi < a < 0.25/mi). 
The dust lies in a spherical shell with a 1/r 2 density distribution and a R out jR[ n = 10 axis ratio between the 
inner and outer radii. While the directly escaping light depends only on the optical depth, thicker shells 
(larger R out /Rin) allow more scattered light to escape than thin shells for the same total optical depth, but 
this is a second order effect and we will hold R ou t/Rin fixed. The model parameters are the luminosity of 
the star, L*, the temperature of the star, T*, the temperature Tj of the dust at R in , and the dust optical depth 
ry at V-band (0.55/im). This is the total (absorption plus scattering) optical depth, r tol = T a j, s + T scat , where 
the scattering optical depth is related to the total optical depth by the albedo, T scat = wr tot . The effective 
absorption optical depth is approximately r e = (t^t^) 1 / 2 = (\-w) l / 2 T tot . At V-band, the standard DUSTY 
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silicate (graphitic) dust models have albedos of wy = 0.86 (0.47), so for the same effective absorption the 
silicate models require significantly higher total optical depths ry. The V-band extinction is approximately 
Ay = 2.5r e .y/lnl0 = 0.41ry (0.79ry) for the silicate (graphitic) model. The radius Rj„ is then a derived 
quantity given the model parameters. We fit the data and estimate uncertainties by embedding DUSTY in 
a Markov Chain Monte Carlo (MCMC) driver. Without mid-IR detections, we cannot determine the dust 
temperature, so we simply considered cases with = 100 K, 500 K, 1000 K and 1500 K. We assume a fixed 
foreground Galactic extinction o f E(B — V) = 0.05 ma g (modeled using an Ry = 3.1 extinction law) and a 



distance modulus of ^ = 30 mag (IFreedman et al.ll2001 



). 



Fraseretal 



While t hey ar e not physically correct models, we do recover the results of I Van Dyk et all (120121) and 
(2012 ) if we mode l our version of the SE P including the mid-I R flux limits but no dust emission 



and using ICardelli et all (|1989l ) extinction laws. Like lVan Dyk et all (|2012t) we find that higher Ry ex t inctio n 
laws provide better fits, but we also find a stellar temperature degeneracy similar to Fraser et all (|2012h . 
Clearly the models have a "degeneracy" direction which is very sensitive to the exact assumptions of the 
model and can lead to eithe r tight or loose limits on the effective temperature. The specific rationale for using 
Ry ~ 4 extinction laws in lVan Dyk et all (|2012r) is not really validjjbut allowing freedom in the extinction 
curve is certainly a more c onservative approach. For example, if we fix Ry = 4.35 and T* = 3600 K to 



Van Dyk et all (120121) , we obtain a lumino sity of \og w L*/L & = 5.15 ±0.05 compared to 5.21 ±0.03. 



match 



If we fix Ry = 3.1 to match 



Fraser et al. 



(120121) . then we find log 10 L»/L Q = 5.10 ±0.05 for T* = 3600 K 
and log lo L*/L0 = 5.43 ±0.05 for = 4400 K, close the values of log 10 L*/L Q = 5.0 and 5.6 they find for 
temperatures of T* = 3550 K and 4450 K using MARCS models with log 10 g = 0. 

Fig.|5]shows two representative DUSTY models of the progenitor for circumstellar silicate and graphitic 
dusts with Td = 1000 K, roughly corresponding to the expected inner edge dust temperatures if the star was 
forming dust at the time of the observations. Formally, the graphitic model is a better fit (x 2 = 5.5 versus 8.8), 



Van Dyk et alj J2012D advocate a high R v based on the finding of [Massev et alj 1(2005) that matching spectroscop ic and B-V 



color e stimates of Ay requires using larger effective value of Rv for the photometric estimate. However, the problemlMassey et al.1 
( 2005) was addressing does not correspond to changing Rv at all wavelengths in a general fit to an SED. iMassev et al.1 J2005I) 
were correcting for the shift in effective filter wavelengths at blue and ultraviolet wavelengths between the O stars usually used to 
build models for the extinction curve and low temperature red supergiants. The effective wavelength of a filter is an average of the 
wavelength weighted by the spectrum over the filter transmission function. Where the spectrum is changing rapidly and cannot be 
modeled as a linear trend, the effective wavelength shifts. This is primarily an issue on the "Wien" tails of stellar spectra, where the 
effective wavelength for a cool star is significantly redder than that of a hot star. For the 3600 K MARCS models, t he B, V and I 



band s hifts relative to an O star are roughly 8%, 3% and 2%, and then become steadily smaller at longer wavelengths. Mass ev et al 



1 2005) derive an effective correction for converting observed B-V colors to V-band extinctions (Av). Mathematically, the observed 
color is (B-V)obs = (R'b-R'v)E(B-V), where R' B and R' v are the Rv = 3.1 extinction law at the shifted, effective wavelengths of the 
B and V bands for the cold star. Thus, getting the correct Av from the observed color corresponds to using R' v = 3A/(R' B -R' V ) so 
that A v = R v (B- V) obs . For the 3600 K MARCS models, the effective B and V-band wavelength shifts then lead to the Ry ~ 4.1 
values reported bv lMassev et al.1 12005). This is a specific correction for the distortion of the B-V colors by the wavelength shifts 
and not a change in the overall extinction curve. At redder wavelengths, the wavelength shifts become steadily smaller and the 
differences between the true and effective (band-averaged) extinction curves become smaller and smaller. Since the B and U band 
fluxes only enter as weak upper limits on the progenitor, it is likely more correct to simply ignore the shifts at V band and longer 
wavelengths than to assume that an underlying Rv = 3.1 extinction curve is globally transformed into aRv — 4.35 extinction curve. 
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Fig. 5. — Graphitic (top) and silicate (bottom) fits to the SED of the progenitor for = 1000 K at the inner 
edge of the dust distribution. The heavy black curves show the SED model, which is comprised of (from 
left to right) scattered, direct and dust emission, as shown by the dashed red curves. The thin blue curves 
show the model for the unobscured SED of the star. The open symbols show the measured luminosities and 
limits. Note the different balance between absorption and scattering for the two dust compositions and the 
modest contribution of dust emission to the K s band. 
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Fig. 6. — Constraints on the visual extinction Ay (top) and stellar luminosity (bottom) as a function 
of stellar temperature T*. The solid contours show the la constraints on two parameters (A% =2.3) for 
the graphitic (left) and silicate (right) models and an inner edge dust temperature of Tj = 1000 K roughly 
corresponding to havin g a dust forming wind at th e time of the observ ations. The points and dashed contours 
show the results from IVan Dyk et all (120121) and Fraser et al.1 (120121) . respectively. The DUSTY model ry 
were converted to A v = 2.5(1 -wv) 1 / 2 ry/ln 10. 
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but the models have essentially identical stellar luminosities, L* = 1O 4 ' 9 L , and temperatures, ~ 3550 K. 
Both models marginally satisfy the mid-IR flux limits at this luminosity, with the graphitic models primarily 
limited by the non-detections at 3.6 and 4.5/mi, and the silicate models limited by the contribution of silicate 
emission peak at 8.0//m. 

Fig. [5] illustrates all three of the basic points about circumstellar dust as compared to interstellar dust. 
First, the optical flux is completely dominated by the scattered emission that is not included in interstellar 
extinction laws. Second, dust emission is quantitatively important to the K s band flux. Because the star is 
cold, dust emission does not dominate the near-IR flux as it would for a hot star of the same luminosity 
and degree of obscuration, but it does partly compensate for the absorption. Third, there are quantitative 
differences between the two dust types in their balance between absorption and scattering and the nature of 
the mid-IR emission. The total optical depths of the models are quite different, ry = 5.9 for silicates and 
Ty = 2.6 for graphite, but the effective absorption optical depths of the two models are very similar, with 
r e V = (\-wy) l l 2 Ty = 2.2 for silicates and 1.9 for graphite. 

Fig. [6] shows the allowed parameter ranges for the stellar lumi nosity an d visua l exti nction Ay as a 



functi on of stellar temperature as compared to the results from Fraser et all (120120 and IVan Dyk et al. 



(120121) . We show the results for Tj = 1000 K, which roughly corresponds to the inner edge dust temper- 
ature if there was a dust forming wind at the time of the observations. Table |2] gives the results for the 
other dust temperatures. They are generally similar except for the Tj = 100 K grap hitic model. The key 
issue for the progenitor is that the preferred solutions of both lFraser et al.1 (|2012l) and Van Dv k et all (120121) 
significantly overestimate the luminosity of the progenitor. Our results agree with Van Dyk et al. ( 2012b on 
the temperature, but the luminosity L* is a fac tor of two lower. T he mid-IR limits essentially preclude the 



higher luminosity and temperature solutions of 



Fraser et al. 



(2012) unless the dust temperature is made sig- 



nificantly lower than Td = 500 K and the dust emission moves out of the IRAC band passes. Even there, only 
the cold Td = 100 K graphitic model allows a significantly higher luminosity. As expected, a large part of 
the difference is that the models based on interstellar extinction are significantly overestimating the overall 
absorption. 



3. Discussion 



Fraser et all (120121) and I Van Dyk et all (12012 ) argue that the progenitor was probabl y the most massive 



yet found for a Type IIP SN, probably at or above the upper limit of (16.5 ± 1.5)M ISmartt et all (120091) 
found in their statistical analysis of the masses of Type IIP progenitors. G i ven t he amount of circumstellar 
extinction, this seemed to match the suggestion by IWalmswell & Eldridgd (120121) that the upper mass limit 
could be biased by increasing levels of dust formation for the more massive red supergiants due to the 
increase in mass loss rates with luminosity. Here we argue that many of these inferences are biased by 
incorrectly modeling circumstellar dust with an interstellar extinction law, thereby overestimating both the 
amount of extinction and the luminosity of the star. When we model the SED using circumstellar dust 
models, the luminosity of the star is L* < 1O 5 L and the mass is M* < 15M , where the downwards shifts 
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Fig. 7. — Extinction A\ as a function of wavelength for silicate (squares) and graphitic (triangles) circum- 
stellar dust with R out /R in = 10 and optical depths of t v = 0.1 (bottom), 0.3, 1.0, 3.0, 10.0 and 20.0 (top). 
Note how the shapes of the curves depend on both composition and optical depth. The curves through the 
points are the model from Table [3] 
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Fig. 8. — Interpolation function = -2.51og 10 F£ for the fraction of the observed light that is di- 
rect emission as a function of wavelength for silicate (squares) and graphitic (triangles) circumstellar dust 
with Rout/Rin = 10 an d optical depths of ry = 0.1 (bottom), 0.3, 1.0, 3.0, 10.0 and 20.0 (top). The curves 
through the points are the model from Table [4] For an input SED S\, the observed SED is S\\0~ 0AAx of 
which 5 a 10-°- 4Aa F^ = S x lO-° A{Ax+G f S) is directly escaping emission and S A 10^- 4A Hl-^A* f ) is scattered 
emission. As in Fig. |7J note how the shapes of the curves depend on both composition and optical depth. 
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are easily understood from the differing physics of interstellar and circumstellar extinction. The visual 
extinction is overestimated because interstellar extinction neglects the contributions of scattered light, and, 
to a lesser extent because of the low stellar temperature, the near-IR (K v ) ext inction is overesti mated by 
neglecting the emission from hot dust. The absence of a mid-IR source noted by Fraser et all (120121) is a key 
point of evidence, since there should have been a detectable source given the proposed, higher luminosities 
unless the dust is cold and emitting at longer wavelengths than the IRAC bands. In this particular case, the 
dust composition has little effect on the inferred properties of the star. 

Unfortunately, without measuring the mid-IR portion of the SED we cannot determine the dust temper- 
ature Td other than limiting it to be lower than the dust destruction temperature (Tj ~ 1500 K). However, the 
wind properties are a strong function of the dust temperature because at fixed optical depth far more mass 
is required if the material is at larger radii and colder temperatures. Ignoring the minor (10%) collections 
from the finite value of R out /Ri„ = 10, the mass loss rate required to support the optical depth is 



M = 4itR in v w T V K v l ~ 10 5 Ri„i5V w ioTv5K v \oo M e vear 1 



(1) 



where R in = 10 15 /?/„i5 cm, v w = 10v w io km/s, Ty = 5tvs and Ky = 10 2 k v ioo cm 2 /g. Table [2] gives the val- 
ues for the individual models. Supporting Ty ~ 5 when Tj ~ 100 K requires mass loss rates of order 
M ~ 10~ 2 - 5 M Q /year and implies total wind masses of order 3M Q that are i mplausible for a star with an 
initial mass of M* < 15M . For comparison, the empirical approximations of Ide Jager et all (119881) predict 
log 10 (M/Mq /year) ~ -5.8 ±0.5 for a red supergiant wind, roughly matching the rates needed if the wind 
is producing dust at the time of the explosion. 

The density of the wind is also tied to the expected phenomenolo gy of the explosion. If w e assume 
the standard p e oc v~ 12 outer ejecta density profile for red supergiants (IMatzner & McKedll999h then the 
expected shock velocity is 



v, = 8200£ 5 7 ( "M 



9/20^-7/20^-1/10^1/10,-1/10 



10 



wio h 



km/s 



(2) 



(e.g. iChevalier & Franssonll2003h . where the total energy of the supernova is E = 10 51 Z?5i ergs, the ejected 



mass is M e = 10M e ioM Q , M = 10~ 4 M_4M Q /year and t\ is the elapsed time in days. A shock expanding 
through a dense medium generates a luminosity of Lj = (l/2)Mv^/v w , which we report in Table [2] for v s = 
5000 km/s and v w = 10 km/s. The emission from the forward shock is usually too hard to be easily detected, 
but assuming that the reverse shock is cooling and that its softer emissions dominate the observable X-ray 
emissions, the expected X-ray luminosity is 



9Mv 3 

L x ~=^- ~ 1.63 x 10 7 £ 



500v v 



l7F 27/20 M -21/20^7/10 v -7/10 f -3/10 LQ 



51 



10 



(3) 



with a temperature of order 1 keV (e.g. IChevalier & Franssonll2003l) . If we model the X-ray fluxes in TableQ] 



using Eqn.[3]and a range of additional absorption from Nh = 10 20 to 10 23 cm 2 , the mass loss rates for epochs 
1, 3 and 5 (which have the smallest uncertainties) are 



M ~ (10 



-6.4±0.4 )10 - 5 .8±0.3 )10 -4.7±0.3) Vh1oM| 



-3/2 -27/14 
el0 c 51 



M /year, 



(4) 
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respectively, with some evidence that the amount of excess abso rption needed above Galactic is increasing 
with time. Such fluctuations and trends are not unusual (e.g. iDwarkadas & Gruszkdl2012T) . but the X- 
ray emission appears to be broadly consiste nt with the presence of a wi nd with roughly the right density 
to explain the extinction of the progenitor. IStockdale et al.1 (120121 ) and lYadav et al.1 (120121 ) report rising 
~ 20 GHz radio fluxes of 0.160 ±0.025 and 0.315 ±0.018 mJy roughly 7 and 13 days after discovery 
that also argue for a s ignificant wind at the time of the SN. If we model the radio emission following 
Soderberg et al.l (|2005l) assuming the ejecta mass is M e = (15- 1.4)M = 13.6M and E51 = 1, we obtain 
estimates of M ~ lO~ 5O v w ioM0/year, although given only two data points at essentially the same frequency, 
the models are not tightly constrained. The Thompson optical depth of the wind is always negligible, since 
it is a small (<^ 1%) fraction of the dust optical depth, but the cold dust solutions would likely convert the 
SN into a Type Iln because the Ha luminosity from recombination is of order 3OOO/?, n i5L0 and increases 
linearly with the distance to the circumstellar material. Thus, while the data is fragmentary, the simplest 
interpretation appears to be that there was a relatively steady M ~ 10 to 10~ 5 MQ/year wind creating the 
obscuration at the time of the SN. 



Since the primary reason for inappropriately using Galactic extinction laws for circumstellar dust is 
almost certainly their ease of use, we supply in Table [3] equivalently easy to use models for absorption 
by circumstellar dust. The problem is somewhat more complex because the extinction depends on both 
wavelength and optical depth, but the absorption in the DUSTY models from the UV to mid/near-IR (0.3/im 
to 5.0fim) can be well modeled by the functional form 



(5) 



for optical depths up to ry = 20. Here A is the wavelength in microns and ry is the total (absorption plus 
scattering) optical depth in the V band. Table [3] provides these models for R out /Rin = 2 and 10 in a format 
where they can simply be extracted from the electronic paper using a mouse and inserted into most numerical 
environments. The fits reproduce the DUSTY results with rms fractional residuals of 1.4% (1.6%) and 1.8% 
(2.2%) for the graphitic and silicate models and R out jR[ n = 10 (2), although this includes some numerical 
rounding errors in the DUSTY models at low optical depths and longer wavelengths. They can be extrapo- 
lated to longer wavelengths relatively safely, but at these longer wavelengths one would almost always also 
need to include dust emission. They should not be extrapolated to shorter wavelengths or higher optical 
depths. Fig. [7] shows the circumstellar extinction of the DUSTY models used to build these interpolating 
functions and the interpolating functions as extracted from Table |3]for R ou t/Rm = 10. In some circumstances 
it will also be useful to separate the direct and scattered light. The fraction of the observed flux that is direct, 
F£ bs , can be well-modeled by 



Gf = -2.51og 10 , A 



F? bs (.T V ) = T v / \l + \ x r 1 



EE 



</ 2 \-i 

aijTy A J . 



(6) 



These interpolating functions are supplied in Table [4] and the quality of the fits is shown in Fig. [8] For 
a source with intrinsic spectrum S\, the observed spectrum is S\\0~ 0AAx , the directly escaping flux is 
S x lQ-° AA xFf s = S x l0-° A(Ax+G f S) and the escaping but scattered flux is S A 10-°- 4A Hl - rahs 



Ff s ). Obviously, 
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this approach can be generalized to other dust compositions or size distributions. We parametrize the mod- 
els by Ty because ry is closely related to the physical properties of the wind (Eqn. [D, while relating it to a 
color (i.e. E(B-V)) would divorce the model from the underlying physics. There is no comparably simple 
means of treating dust emission because it depends critically on the input spectrum. 

These differences are quantitative rather than qualitati ve - in circumstellar dust models the progenito r 
of S N 2012aw is still fairly heavi ly obscured, as found by Fraser et all (120121) and IVan Dyk et al.1 (120121) . 
and lWalmswell & Eldridgd (120121) are correct that circumstellar dust introduces a bias that must be modeled 
whe n considering the s tatistics of SN progenitors, particularly when modeling non-detections in analyses 
like ISmartt et all (120091) . However, since changes in extinction physics exponentially modify quantities 



like luminosities, these differences between the two dust geometries are quantitatively important. We note, 
however, that SN where circumstellar dust will strongly bias inferences about the progenitor have densities 
of circumstellar material that will lead to X-ray or radio emission, as observed for SN 2012aw, because 
the optical depth is proportional to the wind density. For example, combining Eqns. Q] and [3j we see that 
the expected X-ray luminosity is Lx oc ry. This means that the properties of the explosion can be used to 
constrain biases from dust around the progenitor even if the observations of the progenitor are inadequate to 
constrain the circumstellar extinction. 

There are, of course, additional complexities coming from the geometry of the dust around the star that 
can lead to differences in the effective optical depth along the line of sight for direct emission and averaged 
over a significant fraction of the shell for the scattered emission. Our simple models provide two means of 
approximating some of these effects. First, changes in the shell thickness can be used to adjust the balance 
between scattered and absorbed light. Second, the emergent flux can be modeled as a sum of direct light 
with one optical depth, and scattered light with another in order to model differences between the mean and 
line-of-sight optical depths. Independent of these questions, interstellar and circumstellar extinction are 
quantitatively different, and in this case ignoring the differences leads to a significant overestimate of the 
progenitor luminosity and mass. 
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Table 1. Swift X-ray Data 



Obsid 


(Date) 


T ap (sec) 


0.2-10 keV 


0.2-0.5 keV 


0.5-2 keV 


2-10 keV 


Range 


00032319001 


March 21.0 


9151 


1.781 ±0.536 


0.174 ±0.197 


1.445 ±0.458 


0.162 ±0.197 


03/21 


00032319002 


March 23.0 


5856 


0.449 ±0.376 


-0.090 ±0.028 


0.369 ±0.306 


0.171 ±0.217 


03/23 


00032319003 


March 25.0 


8843 


1.528 ±0.497 


-0.053 ±0.018 


1.520 ±0.475 


0.061 ±0.145 


03/25 


00032319004-006 


March 30.4 


22714 


0.618 ±0.229 


0.037 ±0.1 12 


0.323 ±0.146 


0.261 ±0.136 


03/28-04/01 


00032315015-030 


April 12.5 


48188 


0.358±0.136 


0.1 17 ±0.089 


0.059 ±0.058 


0.1 84 ±0.085 


04/01-04/22 


00032315031-042 


June 14.9 


17430 


-0.123 ±0.198 


-0.302 ±0.080 


0.280 ±0.164 


-0.100 ±0.076 


04/26-07/14 



Note. — (Date) is the exposure time weighted average date, where Range gives the range of dates spanned by the epochs. The counts are in units 
of counts/ksec. 



Table 2. Model Fits 



Type 


7rf(K) 


log lo (L*/L ) 


T t (K) 




log| (R,„/cm) 


log 10 M 




log lo (L s /L ) 


sil 


100 


4.99 (4.88/5.12) 


3754 (3483/4168) 


0.87 (0.72/0.97) 


16.97 (16.84/17.10) 


-2.85 (-3.14/- 


-2.63) 


9.16 (8.87/9.38) 


sil 


500 


4.93 (4.80/5.06) 


3600 (3357/3991) 


0.80 (0.48/0.93) 


15.31 (15.13/15.46) 


-4.59 (-5.09/- 


-4.31) 


7.43 (6.92/7.70) 


sil 


1000 


4.91 (4.81/5.01) 


3603 (3385/3976) 


0.79 (0.55/0.91) 


14.80 (14.67/14.93) 


-5.11 (-5.48/- 


-4.86) 


6.90 (6.53/7.15) 


sil 


1500 


4.90 (4.81/4.97) 


3823 (3473/4264) 


0.85 (0.67/0.94) 


14.50 (14.36/14.60) 


-5.35 (-5.67/- 


-5.17) 


6.66 (6.35/6.85) 


gra 


100 


5.14 (4.96/5.30) 


3897 (3528/4308) 


0.64 (0.45/0.74) 


17.40 (17.25/17.53) 


-2.65 (-3.00/- 


-2.43) 


9.36 (9.01/9.58) 


gra 


500 


4.95 (4.82/5.08) 


3528 (3360/3821) 


0.44 (0.13/0.60) 


15.81 (15.70/15.93) 


-4.46 (-4.87/- 


-4.18) 


7.56 (7.14/7.84) 


gra 


1000 


4.89 (4.79/4.97) 


3562 (3372/3910) 


0.43 (0.17/0.57) 


14.97 (14.88/15.07) 


-5.30 (-5.65/- 


-5.07) 


6.71 (6.36/6.95) 


gra 


1500 


4.84 (4.76/4.93) 


3694 (3392/4222) 


0.48 (0.16/0.61) 


14.49 (14.37/14.59) 


-5.74 (-6.17/- 


-5.50) 


6.28 (5.84/6.51) 



Note. — These are the median and 90% confidence intervals of the MCMC results. The implied mass loss rate M is in units of v w \qk v \ oq Mq /year, 
and the shock luminosity L s is in units of (v ! /5OOOkm/s) 3 K l ; 1 |oo L0. 
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Table 3. Graphitic and Silicate Circumstellar Extinction Laws 

Models with R out /Ri„ = 2 

tgral = (0.500446 + 1 .795729 * t- 1 .877658 * t * t + 0.852820 * t * *3 - 0. 141635 * t * *4) 
tgra2= (4.318269- 14.236698 * f+ 13.8041 10 * f* f-5.991369 * f* *3 + 0.959539 *f* *4)// 
tgrai = (-5.114167 + 32.462564*^-26.895305 *t*t+ 10.197398****3- 1.414338 *J**4)//**2 
tgra4 = (3.384105-21. 107633 *t+ 12.229167 * f * f-2.172318 * f * *3-0.149866 *f * *4)// * *3 
t gra5 = (-1 .059677 + 5 .553703 * t- 1 .527415 * t * t - 0.88 1450 * t * *3 + 0.39 1763 * t * *4)// * *4 
tgra6 = (0. 12 1772 - 0.5 1 8968 * t- 0.048566 * t * t + 0.248002 * t * *3 - 0.077054 * t * *4)/Z * *5 

agraphite2 = (tgral+tgra2+tgra3+tgra4+tgrci5 + tgra6) *t* /* *(— 1.375229) 

f si/1 = (0.437549 - 0.446323 * f + 0.648423 * r * f - 0.321970 * r * *3 + 0.055555 * t * *4) 
tei/2 = (-0.486741 + 4.034854*f-5.530127*f*f + 2.711095 *f * *3-0.469112*f * *4)// 
tsil3 = (1.166512- 12.015845 *t+ 14.917191 *f*f-7.058630*f**3+1.204954*f**4)//**2 
tei/4 = (-0.655682 + 1 3 .8495 14 * f- 14.342367 * f * f + 6. 165 157 * f * *3 - 0.984406 * t * *4)// * *3 
f.«75 = (0.169689-4.956815 *f+4.137525*f*f- 1.419899 *f**3 + 0.176166*f**4)//**4 
f«76 = (-0.016829 + 0.582619 *f-0.381166*f*f + 0.083593 *f**3-0.002153*f**4)//**5 

asilicate2 = (tsill+tsil2+tsil3+tsil4+tsil5+tsil6) *t*l* *(-0.642318) 

Models with R oM /Ri„ = 10 
fgral = (0.760499 + 0.879164*f-0.350748*f*f-0.039612*/**3 + 0.034161 *f**4) 
?gra2= (4.061343-7. 166933 *f+2.791544*f*f + 0.214647 *f**3-0.233685*f**4)// 
fgra3 = (-5. 133851 + 16. 344656 *r-4.283100*r*r-1.764900*r**3 + 0.780217*r**4)/Z**2 
fgra4= (3.387184- 10.066016 *f- 1.260999 *r*r + 4.103272 *r**3- 1.160204 *r**4)/Z**3 
t gra5 = (-1 .052057 + 2.479576 * t + 1 .6 1 8868 * t * t - 2.030708 * t * *3 + 0.5 16503 * t * *4)// * *4 
fgra6 = (0.120327-0.214118 * f-0.293914 * f * f + 0.295687 *f * *3-0.071840*f * *4)//* *5 

agraphitel0 = (tgral+tgra2+tgra3+tgra4+tgra5 + tgra6)*t*l * *(— 1 .475236) 

f si/ 1 = (0. 197398 - 0.293417 * t+ 0. 192686 * f * r - 0.041 375 * t * *3 + 0.000902 * t * *4) 
rji/2 = (0.093593 + 2.49 1030 * t- 1 .453387 * f * t + 0.239280 * t * *3 + 0.01 3273 *t**4)/l 
tsil3 = (0.357331 -6.883382 *f+3.239407*f*f-0.198182*f**3-0.121949*f**4)//**2 
tsilA = (0.022567 + 7. 1692 14 * t- 1 .884278 * t * t - 0.728088 * t * *3 + 0.309664 * t * *4)// * *3 
fsi/5 = (-0.065599 - 2. 173 1 30 * f- 0.305525 * f * f + 0.83929 1 * t * *3 - 0.223389 * t * *4)// * *4 
to76 = (0.012188 + 0.216324*f+0.133118*f*f-0.153598*f**3 + 0.036469*f**4)//**5 

asilicatelO = {tsil\+tsil2+tsil3+tsil4+tsil5+tsil6)*t*l* *(-0. 323043) 



Note. — These expressions are designed to simply be grabbed with a mouse from the electronic paper 
and inserted into most numerical environments. The input quantities are t = iv/10 and / = A in microns, and 
the output quantity is A\(ry). They are valid for A > 0.3^m and t v < 20 and should not be extrapolated 
outside this range. No emission by the dust is included in the models. 



Table 4. Direct Emission Fractions for Graphitic and Silicate Circumstellar Extinction Laws 



Models with R out /R in = 1 



ggral = (0.091091-2.031641 *5 + 5.049221 * s* 5-6.399577 * s* *3 + 2.086527 *5* *4) 
#^2 = (-0.705931 + 3.566326*5-40.821335*5*5+48.029644*5**3- 14.815801*5**4)// 
ggrdi = (1.604726-22.761479*5 + 90.223336*5*5-88. 179577*5* *3 + 25.041643*5**4)//**2 
ggraA = (-1 .240894 + 20.383229 * 5-63.820599 * s * 5 + 57.434801 * 5 * *3- 15.475733 * 5 * *4)/Z * *3 
ggra5 = (0.392600-6.8561 17*5+ 19. 106438 * s * s- 16.340312 * s * *3 + 4.240107 * s * *4)// * *4 
ggrab = (-0.044443 + 0.800158 * 5-2.078890 * s * 5+ 1 .716855 * s * *3 - 0.433098 * s * *4)// * *5 

ggraphitel = -(ggral + ggral + ggrai + ggra4+ ggra5 + ggrad) * 5 /(l + / * *(3. 608885)) 
£5«71 =(0.016559-9.435716*5+1.279347*5*5-0.617515*5**3 + 0.351253*5**4) 
g5i72 = (0.054245 + 15.096287 * 5 + 1 .441236 * 5 * 5-0.70455 1 * s * *3 -0.049853 * 5 * *4) // 
g5i73 = (-0.308782 - 6.642689 * 5 - 15 .320068 * 5 * 5 + 12.697774 * 5 * *3 - 3 .365767 * 5 * *4)/l * *2 
gsM = (0.367973-6.190985*5+ 18.372915*5*5- 12.435162*5* *3 + 2.599908 *5* *4)//* *3 
g5;75 = (-0.1 3409 1 + 3 .475205 * 5 - 5 .8 1 0992 * 5 * 5 + 2.902 163 * 5 * *3 - 0.26 1 672 * 5 * *4)/l * *4 
§5//6 = (0.015684-0.468930* 5+0.577899*5* 5-0.158431*5**3-0.045234* 5 **4)/Z**5 

gsilicatel = -(gsill + gsil2 + gsiB + gsil4+ gsilS + gsil6) *s/(l +/* *(4.733179)) 

Models with R ou ,/Ri„ = 10 
ggral = (-0.132213-7.859540 * 5 - 1 .722856 * 5 * 5 + 0.017627 * 5 * *3 + 0.2007 12 * 5 * *4) 
ggral = (0.397806 + 6.85 1899 * 5 + 0.074252 * 5 * 5+ 10.503 128 * 5 * *3 - 3.322806 * 5 * *4)/Z 
ggra3 = (-0.420577-7.736732* 5 + 5.079083*5* 5-1 1.886283 *5**3 + 2.016602*5**4)/Z**2 
ggra4= (0.281626 + 3.202032*5 + 0.592135*5*5+ 1.043499*5**3+ 1.258519*5* *4)/Z**3 
ggra5 = (-0.093566-0.543338* 5- 1.283187*5* 5+ 1.231904*5**3-0.913404* 5 **4)/Z**4 
ggra6 = (0.01 1 334 + 0.026020 * 5 + 0.236427 * 5 * 5 - 0.256076 * 5 * *3 + 0. 140387 * 5 * *4)// * *5 

ggraphitelO = -(ggral + ggral + ggral + ggra4 + ggraS + ggra6) * 5 / '(1 + / * *(4. 39931 1)) 
g«7 1 = (0.041 183 - 10.239763 * 5+3.360815 * 5 * 5-3.994570 * 5 * *3 + 1 .668688 * 5 * *4) 
g5//2 = (-0.127688+ 19.263277*5- 13.077187*5*5+ 19.483457*5**3-8.060155*5**4)// 
g5//3 = (0.141609- 15.254131 *5+ 18.619615*5*5-32.651846*5**3+ 14.198402*5* *4)//**2 
g5//4 = (-0.089791 + 1.766552*5- 13.867558*5*5 + 28.015047*5**3- 12.378082*5* *4)//**3 
gsUS = (0.042445 + 0.565792 * 5 + 5 .72 1465 * 5 * 5 - 10.944883 * 5 * *3 +4.706386 * 5 * *4)/Z * *4 
#5*76 = (-0.006950- 0. 107301 * 5-0.821527 * 5 * 5+ 1 .470719 * 5 * *3 -0.616414 * 5 * *4) // * *5 

gsilicatel0=-(gsill+gsil2+gsil3+gsil4+gsil5+gsil6)* 5/(1 +/**(4 .756651)) 



Note. — These expressions are designed to simply be grabbed with a mouse from the electronic paper and 
inserted into most numerical environments. The input quantities are 5 = (tv/10) 1 / 2 (NOTE THE CHANGE 
FROM t = Ty/10 IN TABLES) and / = A in microns, and the output quantity is G^ s (t v ) from Eqn.0 They 
are valid for A > 0.3/im and Ty < 20 and should not be extrapolated outside this range. No emission by the 
dust is included in the models. 



